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Constraints on the TeV source population and its contribution to 

the galactic diffuse TeV emission 

S. Casanova ^and B. L. Dingus 
Los Alamos National Laboratory, Los Alamos, NM 87545, US 

ABSTRACT 

The detection by the HESS atmospheric Cerenkov telescope of fourteen new 
sources from the Galactic plane makes it possible to estimate the contribution of 
unresolved sources like those detected by HESS to the diffuse Galactic emission 
measured by the Milagro Collaboration. The number-intensity relation and the 
luminosity function for the HESS source population are investigated. By eval- 
uating the contribution of such a source population to the diffuse emission we 
conclude that a significant fraction of the TeV energy emission measured by the 
Milagro experiment could be due to unresolved sources like HESS sources. Pre- 
dictions concerning the number of sources which Veritas, Milagro, and HAWC 
should detect are also given. 

Subject headings: gamma rays: theory 
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1. Introduction 

The Galactic diffuse 7-ray emission is believed to be mostly produced in interactions 
of cosmic rays with the matter and the radiation fields in the Galaxy, the main production 
mechanisms being electron non-thermal Bremsstrahlung, Inverse Compton scatterings off 
the radiation fields and pion decay processes in inelast ic collisions of nucle i and matter. Al- 
though the standard production mechanisms of 7-rays (JBertsch et al.lll993l ) explain generally 
well the spatial and energy distribution of the emission be low 1 GeV, t he model does not 



match EGRET observations of the 7-ray sky above 1 GeV (IHunter et al. 



19971). Manv pos- 



sible e xplanations have b e en pr o posed to account fo r the GeV excess (lAharonian fc Atovan 
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20061 : IStrong I l2006l : IStecker. Hunter fc Kniffen 1120071 ). In any case a major difficulty when 
studying the diffuse emission is to disentangle the truly diffuse emission from that produced 
by unresolved sources. 



The TeV diffus e 7-ray emiss i on fro m the Galactic Plane (jAtkins et al.l 120051 ) and from 
the Cygnus region (JAbdo et al. 1 120071 ) has been me asured by Milagro, a wa ter Cerenkov 
telescope surveying the northern sky at TeV energies (JAtkins et al.ll2000l . l2003l ). The diffuse 



emission measured by Milagro is (7. 3 ± 1.5 ± 2.3) x 10 ^^ photons cm ^ sr ^s ^ for E > 3.5 



TeV in the region 40° < / < 100° (jAtkins et al.ll2005l ). A measurement of the TeV 7-ray 
diffuse flux is important, as it allows us to calculate what fraction of the emission is produced 
through inverse Compton and what fraction is produced through pion-decay mechanism. 
Milagro result, which is consistent with EGRET data above 1 GeV, seems to exclude an 
additional hard spectrum component continuing to above 10 TeV to explain the 60 per cent 
excess of EGRET flux comp ared to vrn- decay production r aecha nism due to the local cosmic 
ray flux (JHunter et al.l 119971 ) . Recently iProdanovic et al.l (120061 ) have argued that no strong 
signal of pion decay is seen in the 7-ray spectrum. The pion decay mechanism would not be 
able to explain the excess above 1 GeV and the Milag ro measurements of the d iffuse flux at 
higher energies reveals a new excess at TeV energies. IProdanovic et al.l (120061 ) investigated 
several possibilities for what might cause the TeV excess, among them possible dark matter 
decay, contribution from unresolved EGRET sources or sources that are only bright in the 
TeV range. An extrapolation of the EGRET sources within the range in Galactic longitude 
of the Milagro observation overestimates the diffuse TeV flux measured. However, there 
could be a population of sources undetectable by EGRET, but contributing to both the 
GeV and TeV excess diffuse emission. Recent observations by the TeV observatory, HESS, 
in fact point to a new class of hard spectrum TeV sources located clos e to the Galactic Plane 



within ±1 degrees of latitude, whose average slope is about E ^-^ (JAharonian et al.l 12005 



2006al ). The HESS detection of high energy 7 rays from fourteen new sources has improved 
significantly the knowledge of both the spatial distribution and the spectra and fluxes of VHE 
7-ray galactic sources. These HESS results make it possible to estimate with unprecedented 
precision the contribution of unresolved sources to the Galactic diffuse emission recently 
extended to TeV energies by Milagro. In fact, the Milagro sensitivity to point sources is 
about 10~^^ photons cm.'"^ s~^ for E > ITeV and it is of the order of the brightest integral 
source flux detected by HESS. Sources like the ones detected by HESS are then unresolved 
by Milagro and contribute to the diffuse emission Milagro measures. 

Here, based on HESS results, knowing the sensitivity and the field of view of an exper- 
iment, we estimate the number of expected sources and their expected VHE 7-ray flux. We 
then evaluate the contribution of unresolved sources to the diffuse 7-ray emission for Milagro. 
The estimate of the diffuse emission arising from unresolved sources will be a lower limit. 
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In fact, HESS has mostly detected extended sources. For sources larger than the angular 
resolution of the experiment the source size determines the sensitivity of the instrument. The 
sensitivity decreases if the extension of the source increases. In Fig.[T]the HESS source fluxes 
and the sensitivity as a function of the source extensions are plotted. A clu stering of sources 



abou t HESS survey sensitivity of 3 percent of the Crab flux reported in (jAharonian et al. 



20051 ) is evident, meaning that there are probably many more to be detected with better 
sensitivity. 
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Fig. 1. — The HESS source f luxes a nd the HESS survey sensitivity versus source extensions 
are plotted lAharonian et al.l (l2006al ). A clustering of the sources about the HESS sensitivity 
limit is evident. 



2. Galactic spatial distributions of PSRs and SNRs and their VHE 

counterparts 

Hereafter the following notations will be used. In terms of the heliocentric distance D, 
the longitude / and the latitude h the galactocentric distance, r, is defined as 



z2 + (r2 + D2_2roL'cos(/)) (1) 

where tq = 8.5 kpc is the distance from the Earth to the Galactic center. The height z over 
the Galactic Plane is 

z = Dsin{b) (2) 

In the inner Galaxy region observ ed by HESS t he existence of 91 SNRs and 389 pulsars 
was already known at lower energies (JGreen I I2004J : Manchester et al.l l2005l ) . Many of these 
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SNRs and pulsars within the inner Galaxy might emit VHE 7 rays but only a few 7 ray 
sources were previously known. Although the HESS Collaboration indicates a clear positional 
coincidence of its sources with a known SNR or pulsar only in a limited number of cases, 
the distribution of Galactic latitude of the seventeen VHE 7 ray s ources detecte d by HESS 
agrees quite well with the distributions of all S NRs catalogued by I Green I ( 120041 ) and of all 
pulsars catalogued by [Manchester et al.l (120051 ) . Assuming, for example, that SNRs are a 
single class of counterparts with isotropic luminosity 1.95 x 10^*^ erg s~^ to the new HESS 
sources and taking for them a simple radiative model, lAharonian et al.l (J2006al ) found that 
the location of these sources favours a scale height of less than 100 pc, consistent with the 
hypothesis that these sources are either SNRs or pulsars in a massive star forming region. 
Therefore although for only a few of HESS sources a firm identification with counterparts at 
other wavelengths exists, there are some suggestions that many of the HESS sources might 
coincide with supernova remnants (SNRs) or pulsar wind nebulae (PWNe). In fact two of the 
HESS sources have SNRs as counterparts, and five of these mo st recently discovered HESS 
sources are associated with pulsar wind nebulae (JFunk 1 120071 ). SNRs are an established 
source class in VHE 7 ray astronomy. Possible correlations between SNRs and unidentified 



EGRET and HESS sources 



(ISturmer &: Dermer 



by (JRomero et al. Ill999l : 



lave been proposed si nce the release of the first EGRET catalogue 



19951: 



Combi et al. 



EsDosito et al. 1119961) . and later for the third EGRET catalogue 



200 ll ). PWNe formed from you ng pulsars with ag e 



less than a million years are considered as potential gamma-ray emitters (JManchester II2005I ). 
Though a young age is not a sufficient condition for a pulsar to generate a PWN. The 
spin-down ener gy loss is the k ey parameter to determine whether a young energetic pulsar 

I2OO3I ). The ratio between 7-ray loud versus 7-ray quiet pulsars is 
Gotthelf I J2003h suggests that all pulsars with dE/dt > dE/dtc = 3.4 x lO^^erg/s 



forms a P WN (IGotthelf 
uncertain. 



are X-ray bright, manifest a distinct pulsar wind nebula (PWN), and are associated with 
a supernova event. By studying t he Chandra data on the 28 most e nergetic pulsars of the 
Parkes Multibeam Pulsar Survey (JManchester II2005I ) iGotthelf I (120041 ) found that 15 pulsars 
with E > 3.4 X lO^^ergs/s are X-ray bright, show a resolved PWN, and are associated with 
evidence of a supernova event. This suggests that about 2.5 per cent of the radio loud pulsar 
have a PWN and might emit 7-rays. 

Supernova remnants and pulsars are the radio counterparts of two of the high energy 
gamma ray candidates, SNRs and PWNe, and their spatial distribution is known from many 



observations at radio wavelengths. The pulsar sur face density apsn 



r), plotted in Fig. El 



is fitted by the following shifted Gamma function (JYusifov fc Kiiciik 

2004J . booeh 



2004; Lorimer et al. 



(7psR{r, z) 



a(-) e' 






(3) 



where a = 41 kpc ^ and b = 1.9 and c = 5.0. For the z distribution we use the exponential 
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function 



npsR{z) = de 



(4) 



where d = 0.75 and e = O.lSkpc (JLorimer et a. 



Fig. El is (JGreenll2004J : ICase fc Bhattacharya 



199 



20061) . The SNR surface density, plotted in 



(ysNR[.r) 



(^osNRsin{^ 



h)e 



-I3r 



ior r < 16. 
for r > 16. 



(5) 



with cTosTVJ? = 1-96 ± 1.38 kpc 



-2 



^2 



17.2 ± 1.9 kpc, 00 = 0.08 ± 0.33 and (3 = 0.13 ± 0.08. 



For the z distribution we use the exponential function 

1^1 



nsNR[z 



de 



(6) 



where d = 0.58 and e = 0.083kpc fIXu et al. 1120051 ). Assuming the SNR and PSR distribu- 
tions plotted in FigJJl the Galaxy should contain 217 SNRs and 13176 PSRs. 
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Pu 



sar a nd SNR s urface d istrib utions versus galactocentric radiu s r from 



Lorimer et al. I (120061 ) and from iGreen I (12004 ) and ICase fc Bhattacharya I (Il998l ). respec- 
tively. 



3. Number-intensity relation for HESS sources 

In order to perform a study of the collective properties of HESS source population the 
number-intensity relation, logN{> S)—logS, is here used. The number-intensity relation has 
the advantage of using the flux data without any assumption on the distance and luminosity. 
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Fig. 3. — Pulsar and SNR surface distributions versus the heliocentric distance D for the 
regions of the Galactic plane surveyed by Milagro ( 40° < / < 100° ) and by HESS ( -30° < 
/ < 30° ). 



In fact for many of HESS sources the location and luminosity are unknown. The number- 
intensity relation also gives information on the geometry of the volume in which the sources 
are contained if a uniform source distribution is assumed. Deviation from a simple power 
law A^ oc 5*^° with slope -1 (for the case of a thin disk) means that the sources or their 
luminosity function are not spatially uniformly distributed. 

The major difficulties in the study of the collective properties of the HESS source popu- 
lation consist of the limited number of sources detected and the relatively small range of flux 
covered by the survey. Also, the HESS survey of the Galaxy was not performed with uniform 
sensitivity. In fact, whereas the sensitivity of the survey of the Galactic Plane in Galactic 
latitude is rather flat in the region between -1.5 and 1.5 degrees, its effective exposure and 
therefore its sensitivity is not uniform in longitude. Longer observation times were dedicated 
by HESS to locations in the Galactic plane close to where three sources, HESS J1747-218 
and HESS J1745-290 (in the Galactic Center), and HESS J1713-397, were already known. 
The average sensitivity of t he survey as a function of the longitude and the latitude are 
shown in Fig. 2 and Fig. 3 of lAharonian et al.l (l2006al ). respectively. In some locations of the 
G alaxy the survey was doii e at peak sensitivity of 2 per cent of the Crab flux. From Fig. 3 
of lAharonian et al.l (l2006al ) one can deduce that in order for our sample to be complete, 
only sources detected with more than 6 per cent the Crab flux within —2° < Z < 2° can 
be included. Table 1 lists the eleven sources included to plot the number-intensity relation. 
HESS source fluxes have small statistical errors, whereas their systematical uncertainties 
are about 30 per cent of the absolute value of the flux. In order to account for the large 
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SOURCE NAME 


FLUX ABOVE 200 GeV 


HESS J16U- 518 


57.8 ±7.7 


HESSJIQIQ - 508 


43.3 ±2.0 


HESS J1632 - 478 


28.7 ±5.3 


HESSJ1634: - 472 


13.4 ±2.6 


HESSJ16A0 - 465 


20.9 ±2.2 


HESS J1702 - 420 


15.9 ±1.8 


HESSJ1804: - 216 


53.2 ±2.0 


HESSJ1813 - 178 


14.2 ±1.1 


HESSJ1825 - 137 


39.4 ±2.2 


HESSJ1834: - 087 


18.7 ±2.0 


HESS J1837 - 069 


30.4 ±1.6 



Table 1: HESS sources included in the logN-logS with fluxes above 6 per cent of the Crab flux, 
All flux es for E > 200 GeV are in units of in 10^^"^ photons cm^'^ sr^^ s~^ lAharonian et al, 
f l2006ah . 



systematic uncertainties in the fluxes a Monte-Carlo based procedure is adopted to generate 
the errors for the logN[> S) — logS diagram, plotted in Fig. HI A random number spanning 
an interval equal to the systematic and statistical errors in quadrature around the flux value 
i s generated. By repeating the procedure many times standard deviations can be evaluated 
( iBignami fc Caraveo Ill980l ). These standard deviations are the error bars in the logN — logS 
diagram plotted in Fig. HI The power law fltting curve has a slope — 1.0 ± 0.1 with reduced 
X^ = 0.5 for the HESS sample for integral fluxes bigger than ^i = 12 x 10^^^ cm^ s~^, which 
corresponds to 6 per cent of the Crab flux. The number-intensity relation is 



C (-l.OitO.l) 

Ar(>S) = (152±41)(-) 

•JO 



(7) 



where Sq = 10 



-12 



cm 



2 „-l 



s . The reduced X for a flt with slope -1.5 would be 3.3. By assuming 



a tridimensional isotropic distribution of sources the slope of the number-intensity relation 
should be -1.5. The fact that the slope is about -1 means that the source decrease in the flux 
with distance is counterbalanced by the density of sources depending on the inverse square 
of the distance. The HESS-like sources are distributed in a thin disk, in a volume that is 



larger than the visibihty hmit. 
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Fig. 4. — The number-intensity relation, logN{> S) —logS, for HESS source population. In 
order to have a complete sample, only sources detected above 6 per cent of the Crab flux 
are included. The population of sources is distributed in a thin disk. 



3.1. Predictions 



The 7-ray flux due to the HESS source population above 6 per cent of the Crab flux in 
the region -30° < / < 30° and -2° < / < 2° is 



,-2 „-l 



F{E > 200GeV, -30 < / < 30, -2 < 6 < 2) = / N{> S) dS = 2.5x10-^° photons cm.-^ s 



where ^2 = 57.8 x 10 



'12 



cm' 



s ^ is the maximum flux detected by HESS from a source. 



From the PSR and SNR distributions plotted in Fig. ([2]) there are 84 SNRs and 5324 
PRSs in the region —30° < / < 30° and —2° < b < 2° , whereas in the Milagro region 
(40° < / < 100° and -5° < 6 < 5°) there are 36 SNRs and 1358 PSRs. Assuming that 
supernova remnants and pulsars are the radio counterparts of high energy gamma ray sources, 
SNRs and PWNe, the Milagro region has 26 percent of the sources which are in the HESS 
region. Also, the Milagro Galactic diffuse emission is measured for a threshold energy of 
about 3.5 TeV, whereas the flux from unresolved sources calculated in Eq.([S]) refers to the 
HESS threshold energy of 200 GeV. In order to estimate the contribution of unresolved 
sources to Milagro diffuse emission the flux in Eq.([8]) has to be corrected for the different 
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threshold energy. All HESS sources are fitted by a power law spectrum 

$(i?,r) = $o(Y^) , (9) 

and the average spectral index is F = 2.32. By correcting the flux in Eq.(IH]) for the Milagro 
threshold energy the flux which HESS-like sources contribute is 

F{E > 3.5TeV, 40 < / < 100, -5 < 6 < 5) = 0.26 F{E > 200 GeV, -30 < / < 30) 

F (-r+i) 

/ -'^Milagro \ 

Ehess 
= 8.3 X lO^"*^^ p/iotons cni~^sr^"'^ s^^ . 

(10) 

The contribution of HESS source population amounts to at least 10 per cent of the diffuse 
flux which Milagro measures above 3.5 TeV. This is a lower limit for the contribution of 
unresolved sources to Milagro diffuse emission, as only sources above 6 percent of the Crab 
flux were taken into account to estimate it. 

The Galactic diffuse emission measured by Milagro can be used to constrain the mini- 
mum flux Smin below which the logN-logS plot becomes flat in order not to overproduce the 
Milagro flux 

/ dS—-<7.3x 10"" photons cm^^ sr"^ s"^ . (11) 

la . dS 

Smin cannot be less than 1 x 10^^^ photons sr~^cm~^ s~^ in order not to violate the constraint 
in Eq.dTTD. 



From the number-intensity relation for HESS source population it is possible to deduce 
the number of sources which HESS, VERITAS, Milagro and HAWC will detect. The number 
of SNRs and PWNs expected for HESS if its entire fleld of view is scanned with a uniform 
sensitivity of 2 per cent of the Crab flux above 200 GeV is about 43 ± 10. If VERITAS will 
survey the Northern sky in the region between 30° < / < 220° and —10° < b < 10° reaching 
the level of 1 per cent of the Crab flux above 100 GeV, it should detect approximately 35 ±10 
sources. Milagro has previously observed the Northern sky with a sensitivity equal to about 
65 per cent the Crab flux at 3.5 TeV. At this level of sensitivity one shoul d expect no de t ection 



of sources for Milagro and indeed no detection of sources was claimed in lAtkins et al.l (120051 ) 
with the data accumulated after the flrst three years of operation. Now after six years of data 
have been accumulated Milagro has reached a sensitivity equal to about 20 per cent of the 
Crab flux at 20 TeV median energy and should be able to detect 3 ± 1 sources. The Milagro 
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CoUaboration has recently published the survey of the Northern sky in the region between 
30° < / < 220° and -10° < 6 < 10° at a thresh old of 20 TeV and a threshold sensitivity of 



about 20 percent of the Crab detecting 4 sources I Ab do et al. I (120071 ). The number of sources 
detected is in agreement with our predictions. After only one year of operation the proposed 
experiment HAWC will have surveyed the region of the sky with longitude 5° < / < 110° 
and 130° < / < 250° and latitude -10° < 6 < 10° at 30mCrab sensitivity above 1 TeV and 
should have detected 19 ± 5 HESS-like SNRs and PWNe. The predictions given here for 
Milagro and HAWC are valid unless the source spectra steepen or cut-off below the Milagro 
threshold. However, for the sources detected by HESS the spectrum is predominantly well 
characterized by a single power law. 



4. Alternative method to estimate the contribution of unresolved sources to 

the diffuse emission 

Thanks to the logN-logS relation we obtained a lower limit for the contribution of 
unresolved HESS-like sources to the diffuse emission measured by Milagro. Here we will 
assume that the density of VHE 7-ray source candidates, SNRs and PWNe, follows the 
volume density of SNRs or of pulsars in the Galactic Plane as observed at radio wavelengths 
and their luminosity function is a power law, and we will then estimate their contribution 
to the Milagro diffuse emission. 

Since the data on HESS sources are too sparse to constrain their luminosity function, 
we will leave it as a parametrised input. The luminosity function $(L) will be a power law 
with different indices a varying between -1 and -2 

^{L) = ^ = c{'^) . (12) 

The assumed luminosity function will then be compared with the HESS source counts to fix 
the normalisation c. In Eg. (1121) L^^ = 1 x IQ^'^erg/ s. The range in luminosities for the HESS 
sources, for which the distance and thus the luminosity is known, varies between IQ^^erg/ s 
and IQ^^erg/ s. In fact, most sources of 7-ray in the Galaxy are located close to the plane 
of the Galaxy, within a region which extends from Dmin = 0.3 kpc up to D^ax = 30 kpc 



(jSwordy II2003I ). The range in luminosity for the HESS sample can then be found from the 
HESS sensitivity (we assume 6 percent of the Crab flux) and the maximum flux detected by 
HESS, which are respectively 



r-1 



^7mm = T^ 7; ^th 4 TT L'mm^ fmin = 3 X lO^^Crg/s 
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(13) 

where Eth is the detector threshold energy and F is the spectral index if the 7-ray emission 
is a power law. By requiring that the following integral over HESS field of view and over the 
luminosity range given in Eq.( fT3l) gives the number of sources HESS detects above 6 percent 
of the Crab flux we will obtain the normalisation factor c 

N(L^ . ,L^ ,V) = dV ^^^^Wes^ 

JVhESS JL^rmn ^^^^ 

dV / dLpsources{r,z)c{—^ ) 

VhESS -^L^rmn ^0 

30 /-so 1-2 fL^^^^ T « 



dD D I dl dbsin{b) psources{l,b, D) / dL c{ , 

0.3 i-30 J -2 JL^rmn ^T 

(14) 
In Eq.l TT^ Psources is the density of SNRs and PWNe. 



4.1. Predictions 

The integral flux above Eq from resolved and unresolved sources from a region in the 
Galaxy which extends from Dmin to Dmax in heliocentric distance, from hmin to hmax in 
latitude and from Imin to Imax in longitude is 

F{E > Eq, Imin < I < ^max, bmin < b < bmax) = / dL^ "^(-^7) -^(-^7, D) dNsources{r, z) 

Jo 

dL-,4(L,)F(L,.Z))/dV'.„„„„(r),w„.(.) 









dL^^{L^)F{L^,D) / dDD'^ 

J Dmin 

/ dl dbsin{b) p{l,b,d) 

(15) 
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where 



If the luminosity function defined in Eq. flT2|) has slope a = —1 the total flux from 
resolved and unresolved PWNe divided over the HESS field of view if 2.5 per cent of radio 
loud pulsars emit 7-rays is 

Fpwn{E > 200GeV, -30 < / < 30, -3 < 6 < 3) = 1.3 x 10"^ p/iotons cm'^ sr^^ s~^ . (17) 

The total flux from the Galactic SNR population is 

FsNRiE > 200GeV, -30 < / < 30, -3 < 6 < 3) = 1.1 x 10'^ photons cm~^ s'^ st'^ , (18) 

For a = —1 the integral flux expected for the Milagro field of view from HESS-like 
sources assuming a pulsar shaped surface density if 2.5 per cent of radio loud pulsars emit 

7-rays is 

Fpwn{E > 200GeV, 40 < / < 100, -5 < 6 < 5) = 4.6 x 10-^>/iotons cm^^ s"^ sr"^ . (19) 

while the contribution from SNRs shaped surface density is 

FsNuiE > 200GeV, 40 < / < 100 - 5 < 6 < 5) = 4.4 x IQ-^^ photons cm'^ s"^ sr"^ . (20) 

The fluxes in Eq. (1191) and (I2U]) are lower than the HESS fluxes from resolved and unresolved 
sources. In fact, as shown in Fig|3l HESS observes the inner region of the Galaxy, whereas 
Milagro field of view is more spread toward the outer regions in the Galaxy, where the 
number of sources is substantially lower, so a lower contribution from unresolved sources is 
expected for Milagro's region of the Galactic plane than for HESS. 

The integral flux from unresolved pulsars if 2.5 per cent of radio loud pulsars emit 7-rays 
corrected for the Milagro threshold becomes 

Fpwn{E > 3.5 TeV, 40 < / < 100) = Fpwn{E > 200 GeV, 40 < / < 100) 

F (-r+i) 



X 



Ehess 
= 3.6 X 10^^^ photons cra~^si^^ s^^ 



(21) 
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whereas the contribution from SNRs is 

FsNRiE > 3.5 Tel/, 40 < / < 100) = Fpuisar{E > 200 GeV, 40 < / < 100) 

F (-r+i) 

Ehess 
= 3.4 X 10^^-^ photons cra~'^ST^^ s'^^ 



(22) 



where EMUagro and E^ess are the Milagro and HESS energy thresholds, respectively. Assum- 
ing a slope a = —1. for the luminosity functiony and summing the two contributions from 
HESS-like SNRs and PWNe, the emission due to unresolved sources to the Milagro diffuse 
emission is comparable to the diffuse emission itself. If the slope of the luminosity function 
a = —1.5 the contribution of unresolved SNRs and PWNe to the diffuse emission measured 
by Milagro is about 5 percent and for a = —2 this contribution becomes negligible, which is 
in disagreement with the lower limit of 10 percent for the contribution of HESS-like sources 
to the VHE diffuse emission previously found. Thus the slope of the luminosity function for 
HESS-like sources is constrained to be —1 > a > —1.5. 



5. Conclusions 

The number-intensity relation and the luminosity function for the HESS source popula- 
tion were investigated using the assumption that HESS sources are distributed as PSRs and 
SNRs detected at radio wavelengths. In order for the chosen sample of sources to be complete 
only the HESS sources with fluxes above 6 percent of the Crab flux were taken into account 
to derive the number-intensity relation. The contribution of unresolved HESS-like sources to 
the diffuse emission measured by Milagro was also estimated. Using the logN-logS relation 
for the HESS sample of Galactic 7-ray emitters at least 10 per cent of the diffuse emission 
at TeV energies is estimated to be due to the contribution of unresolved HESS-like sources. 
This result is a lower limit for such a contribution because we have taken into account only 
sources detected above 6 per cent of the Crab flux and because HESS sensitivity gets worse 
for extended sources, meaning that some extended sources might have been missed by HESS. 
Using the logN-logS relation we have also predicted the number of HESS-like sources which 
VERITAS, HESS and HAWC should detect during their survey of the sky. An alternative 
procedure to evaluate the contribution of unresolved HESS-like sources to Milagro diffuse 
emission gives the diffuse flux due to unresolved sources comparable to the diffuse emission 
itself. We flnally constrained the slope of the luminosity function. The main uncertainty of 



-14- 



this calculation consists in assuming that the distribution of 7-ray sources follows the dis- 
tribution of either pulsars or SNRs observed in the radio. In particular, in order to predict 
how many PSRs observed in the radio have a PWN and are possible gamma ray emitters we 
used the result that the spin-down energy loss dE/dt > dE/dtc = 4 x lO^^erg/s for a young 
energetic pulsar to form a PWN. In this respect we have ignored the existence of pulsars, 
such as Geminga, which are 7-ray loud, yet not observed in the radio. 

New observational results support the hypothesis that a population of unresolved sources 
contribute significantly to the emission at very high energy. Milagro has recently reported the 
discovery of TeV gamma ray emission from the Cygnus Region of t he Galaxy, which ex ceeds 



the predictions of conventional models of gamma -ray production (jAbdo et al. 1120071 ) from 



the s ame region in the Galax y where the Tibet Array has detected an excess of cosmic 



rays (jAmenomori et al. II2006I ). Thanks to its improved sensitivity Milagro has also better 



i maged the whole Northern sky and discovered four sources and four source candidates 



( lAbdo et al. 1120071 ). HESS has seen very high ener gy emission spatially co rrelated with 



giant molecular clouds located in the Galactic Center (jAharonian et al.ll2006bl ). The energy 
spectrum measured by HESS close to the Galactic Center is E~'^'^, significantly harder than 
the ii^~^-^ spectrum of the diffuse emission and equal to the average spectrum of the HESS 
source population. The emission from the Galactic Center might possibly unveil a cosmic 
ray accelerator. 

To draw more definitive conclusions about the very high energy 7-ray sky, new obser- 
vations are of fundamental importance. New hints will be provide by both MAGIC and 
VERITAS, which already survey the Cygnus Region. Finally GLAST will investigate the 
window of energy between 10 MeV to 300 GeV, covering the energy gap left between EGRET 
and the ground-based low threshold gamma-ray observatories. 
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